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A series of lenses whose focusing properties are alternately convergent and 
divergent (alternating-gradient focusing) is of potential interest in the guid- 
ance of light waves, and has previously been used to focus electron beams and 
high-energy particle streams. New information is provided herein on such 
focusing for the case of equal focal length f (but alternating-gradient) lenses 
equally spaced a distance L. 

The alternating-gradient system formed by adding diverging lenses between 
the lenses of an all-converging sequence of lenses is found to have the same 
stability condition as the original system for < L/f < 2. A physical argu- 
ment leads to the conclusion that weaker divergent lenses would also leave the 
stability criterion unchanged. 

The focusing effect of the alternating-gradient system is surprisingly close 
to that of an all-converging lens system. After the focal length of each has been 
adjusted to an optimum value, the ray departure from the system axis is only 
1 .07 times as great for the alternating-gradient system as for an all-convergent 
lens system with the same spacing of convergent lenses. 

For weak lenses (i.e. 2f/L » /) the output ray departure due to input ray 
displacement is independent of both the focal length and spacing of the lenses, 
and is independent of lens spacing but proportional to focal length for input 
ray slope. 

Both the alternating -gradient system and all-convergent lens focusing ar- 
rangements exhibit discontinuities in the maximum ray displacement ver- 
sus focal strength relation. 

Viewed over-all, alternating-gradient focusing for light guidance does a 
surprisingly efficient job and may be advantageous over all-convergent lens 
systems if the alternating-gradient arrangement has structural or economic 
advantages. 
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I. INTRODUCTION 

In research on guidance of light waves for communication we are con- 
sidering use of a sequence of lenses of alternately convergent and diver- 
gent types. For example, a guidance system using tubular thermal gas 
lenses might employ continuous flow of gas through a tube whose walls 
are alternately warmer and cooler than the gas within. Thus, the mech- 
anism used to cause the focusing may have the alternating character, 
and the question comes to the fore — how well can one focus with such 
a structure as compared to the use of a sequence of all -convergent lenses? 

Alternating-gradient focusing has previously been used on electron 
beams 1 and on particle accelerators. 2 The present study discovered an 
error in the previous determination of stability conditions and revealed 
some little known but interesting properties of alternating-gradient focus- 
ing. A comparison is made with conventional focusing using all-conver- 
gent lenses. 

II. ANALYSIS OF ALTERNATING-GRADIENT SYSTEMS 

One might wonder whether a series of equal-focal-length and alter- 
nately converging and diverging lenses would give any net focusing at 
all, since the average dielectric constant along all paths parallel to the 
axis would be the same. It is well known, however, that a divergent 
lens followed by a convergent lens of equal focal length, spaced a finite 
distance less than the focal length, gives a net converging lens, and the 
same is true if the order of the lenses is reversed. Thus, a net focusing 
is to be expected for an infinite series of such lens pairs. 1 

We consider a sequence of alternating convergent and divergent lenses 
equally spaced a distance L and of equal focal lengths, /. We follow the 
method of analysis used by Pierce. 1 There are two cases to cover, one in 
which the first lens of the array is a divergent lens (obtained by starting 
at n = in Fig. 1) and the other in which the first lens of the array is a 
convergent lens (obtained by starting at N = in Fig. 4, below). 

With reference to Fig. 2, and taking the input ray at plane a to have a 
slope r a ' and a displacement r a from the longitudinal axis, the output ray 
from the lens at plane b will be 

r b = r a + Lr a ' (1) 

r b ' = r a '+ (\/})r b . (2) 

At plane c this ray will be described by 

r c = n + Lr b ' (3) 

(4) 



A LTERNATING-GRADIENT FOCUSING 



1743 



L__ 



, L 



t 1_ 



« L ^ 



, 1_ , 



Fig. 1 Sequence of lenses — first case. 



Eliminating n, and r h ' from (3) and (4) gives 

L 



Hence, vvitli reference to Fig. 1 , we can write 

r„ + , = (l + j) r„ + L ('2 4 ': ) ,-,' 



r»+i + r„ = 0. 



These two equations lead to 

*•-[-©: 

The solution to (9) is 

/■„ = -1 cos nd + # sin 7?2 



(5) 
(6) 

(7) 
(8) 

(9) 

(10) 



.i- L ... 



Fig. 2 — Lens subsection for Fig. 1. 
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where 

— cos 



•u-\m\ 



and where A and B are constants to be determined. 

Equation (11) differs from the corresponding equation on page 200 of 
Ref. (1), which is believed to be in error. The correct condition for 
stability, from (11), is 

< \{L/ff < 2 (12) 

or 

< L/f < 2. 

We put in the boundary conditions, at n = 

■/■„' = r„' (13) 

r B = r . (14) 

We make use of a general theorem* stating the orthogonality of the effects 
of r and n' and seek a solution with those quantities as factors. This 
leads to the following form for r„ , using (10), (7), (13) and (14) : 

r n = r fa cos (nd - ipi) + r 'L fa sin nd (15) 

where 



*-b= 



W2f)_ 

<pi = | cos 



(10) 



-HrM d7) 



, _ [(2//L) + 1] _ [2 + (L/f)) ( . 

2 ~[l-\(L/m sin0 "• 

The general form of fa and k 2 versus L/f is shown in Fig. 3. Further dis- 
cussion will be postponed to a later point in this paper. 

We are also interested in the displacement r m at the output of the 
?nth diverging lens (Fig. 1). Using the relation 

r m = r n + Lr n ' (19) 

and using (7) for r„', (15) for r n with appropriate trigonometric rela- 
tions, it can be shown that 

r m = r k 3 cos (md - <p 3 ) + n'L fc 4 cos (md - <p 4 ) (20) 

* See Appendix B. 
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Fig. 3 — Coefficients relating input ray slope and displacement to ray displace- 
ment at the nth, with, Nth, and 1/th lenses of the alternating-gradient lens sys- 
tems of Figs. 1 and 4. 



where 



|_1 + (L/2/)J 



+ (WW. 
k 4 = 2f/L 



(21) 



I — 1 7 — 1 

ip-i = cos Its 



(22) 
(23) 

(fi = I COS K4 (24) 

and is again defined by (11). Plots of k 3 and k A are given in Fig. 3. 



I -1 7 -1 

= I COS IC4 
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Equations (15) and (20) give the ray displacements at any lens in the 
system when the input is at a plane adjacent to a converging lens (i.e. 
at n = 0, Fig. 1 ) . Before discussing interesting features of such ray prop- 
agation we will give the corresponding solutions for the case where the 
input is adjacent to a diverging lens (i.e. at TV = 0, Fig. 4). 

With reference to Figs. 4 and 5, and following a derivation similar to 
that carried out in connection with equations (1) to (9), it is found that 

r N+i - [2 - (L/f) 2 } r N+1 + r N = 0. (25) 

Note that (25) is identical to (9) and with the change of n into N the 
solution for (25) is again (10) and (11). When the initial conditions are 
put in, &t N = 0, r N = r Q and r N ' = r ', we get 



where 



r N = r k r cos (Nd + <p b ) + r 'Lk 6 sin Nd (26) 

= (L//)[2 - (L/0 m I" 2 "I* 

sine Ll + (L/2/)J U7; 

<p 5 — | cos -1 kf 1 | (28) 



K2//L) - 1] [2 - (L/f)} 
K6 [l--j(L/OT" sine ' 



(29) 



Note that A- 5 is identical to k* in (21). 
For r M we find 

r M = r k 7 cos (Md + ^7) + r 'Lk s sin (Md - <p s ) (30) 

where 

*• - [i^iml (31) 

<p 7 = I cos -1 k 7 ~ x I (32) 

A-8 = 2f/L (33) 

<p 8 = I cos -1 /c 8 _1 I . (34) 

Note that A' 8 and <p s are identical to k t and ip* , and that A: 7 is identical to 
ki . Plots of h , kt , ky , k s are given in Fig. 3. 

III. RELATIONS FOR A SEQUENCE OF CONVERGING LENSES 

For comparison purposes we will want to refer to the case of a se- 
quence of identical convergent lenses equally spaced. The analysis is 
similar to that above for the alternating gradient lenses. The results are 
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Fig. 4 — Sequence of lenses — second case. 

as follows. Let the lens spacing he s, and the focal length be/. Then the 
displacement r p at the pth lens is 



where 



r P = r k 9 cos (p8 — <&) -f r 's k l0 sin p8 

a- = r 4//s t 

_(-if/s) - lj 

s _1 kf 



Va = I cos Afe ! 



A' 1 1) — 



1 



11 " («/4/)]* sin 8 
8 = cos - ' [1 - (s/2/)] 
s// = 2(1 - cos 8). 



(35) 

(36) 
(37) 
(38) 

(39) 

(40) 



The system is stable in the sense that an input displacement r or slope 
/V will remain bounded as p is increased if 

< s/f < 4. (41) 

Fig. 6 shows the values of /,- 9 and /wo for comparison to Fig. 3. 

ABC 



— L ■ 4* L 

Fig. 5 — Lens subsection for Fig. 4. 
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Fig. 6 — Coefficients relating ray slope and displacement to ray displacement 
at the pth lens of an all-convergent lens system;/ = focal length, s = lens spacing. 



IV. STABILITY COMPARISON 

Since the converging lenses in Fig. 1 or Fig. 4 are spaced a distance 2L, 
a comparison of (41) and (12) shows that the alternating-gradient system 
formed by adding a divergent lens halfway between the convergent lenses 
of an all-convergent lens system has the same stability condition as the 
original system. At a later point it will be shown that this is reasonable 
physically. 

V. THE WEAK LENS CASE 

When the lenses are weak, i.e., when 2f/L » 1, the general expres- 
sions may be simplified to show some remarkable properties of alternat- 
ing-gradient focusing. When the first lens is a diverging one, (15) and 
(20) yield 

r m = r n = r \/2 cos [nd - (x/4)] + r '2/sin nd (42) 

and when the first lens is a converging one, (2(5) and (30) yield 

,. v _ rM = ro V2 cos [N0 + (t/4)] + r„'2/sin Nd. (43) 

These expressions show that for an input ray displacement without 
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slope the maximum displacement of the transmitted ray (as n or N 
varies) is \/2 times the input ray displacement, independent of both 
focal length and lens spacing! Also, for an input ray of zero displacement 
but finite slope r ', the maximum displacement of the transmitted ray 
is 2/V'o', independent of lens spacing L. The angle 6 is dependent on / 
and L and goes to zero as / — > oo . 

In an all-convergent lens system the similar condition 4//s » 1 leads 
to [from (35)] 

i-p = r cos p8 + r ' y/fs sin p8. (44) 

In comparing the alternating-gradient system to the all-convergent 
lens system for weak lenses, we see that for an input ray with zero slope 
/•(/ but finite displacement, ? - o , the maximum output displacement for 
the alternating-gradient system is y/2 times that of the all-convergent 
lens system. For input ray displacement r = but finite rj, we see that 
the maximum output displacement for the alternating-gradient system 
is larger than for the all-converging lens system by the factor [see (42) 
and (44)]: 

Jfe " (4//s) ' (45) 

Our assumption of weak lenses made 4f/s ^> 1, so (45) is a factor of two 
or more. 

In this weak lens case both 6 and 8 are small angles, and from (11) 
and (39) 

9 ^ L/f (46) 

»«(•//)*. (47) 

Using the case of s = 2L, which is the alternating-gradient system formed 
by adding a diverging lens between the lenses of an all-convergent lens 
system 

6/8 = (s/4f)\ (48.) 

Since 4f/s » 1 by our weak lens definition, 6/8 is less than unity and the 
period of the alternating-gradient system encompasses a great many 
more convergent lenses than does the all-convergent lens system with 
the same spacing of convergent lenses. This is as would be expected. 

VI. OPTIMUM FOCAL LENGTHS 

We now inquire as to whether there is a best value for the lens strength 
in order to minimize output ray displacement. On the assumption that 
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the sine and cosine terms of (15), (20), (26), (30) and (35) go through 
unity for some number of lenses, the question is whether or not the 
coefficients fci , /.- 2 • ■ • A- 10 have any minima. 

For the all-convergent lens system Fig. 6 illustrates that fc 9 has no 
useful minimum but that fe u , relating input ray slope to output ray 
displacement as in (35), does have a minimum. By setting 



% (ho) = 
df 



we find 



/ 



= 2 



(49) 



(50) 



opt 



at which condition fc 10 = 1.0. We note that the displacement r p due to 
ro' is n/skio sin p8, so we have a minimum in this displacement when fc w 
is a minimum provided sin p8 goes through unity for some number of 
lenses p. This is the most typical case, but there are notable exceptions. 
Suppose, for example, that 5 of (39) is tt/3, corresponding to s/f = 1; 
then p8 = tt/3, 2t/3, t, 4x/3, etc., as illustrated in Fig. 7, and | sin p8 | 
never exceeds sin 5. Hence the maximum value of r 's k i0 sin p8 is r 's 
for s/f = 1. It is shown in Appendix A that there is an infinite series of 
such discrete values, but the largest departure of the maximum value of 
rn's/cio sin p8 from k lQ is 15 per cent, occurring at s/f values of 1 and 3, 
as illustrated in Fig. 14 (see Appendix A). 

Turning now to the alternating-gradient system, the only coefficient 
having a useful minimum is lc 2 of (15), relating input ray slope to ray 
displacement at the converging lenses of Fig. 1. We find the minimum in 
fc 2 by setting 

' l (Jfc 2 ) = (51) 



df 



which leads to the equation 

(L//) 8 + 4(L//) 2 -8 = 0. 



(52) 
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Fig. 7 — Diagram of sin pS, p = 1,2,3 ■ • • , when 5 = tt/3. 
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The appropriate root of this equation is 

L 
f 



= 1.237 (53) 

opt 



at which we calculate 6 = 76.4°, and /c 2 = 3.33. Fig. 3 shows this mini- 
mum is sharper than the corresponding one for few of the all-converging 
lens system, Fig. (5. We note that fcj sin nO of (15) contains the sin nd/s'm 
6 factor, so once again (as described in Appendix A) for 6 = tt/3 and other 
values, the maximum value of /i 2 sin nd will be somewhat less than the 
value of ko . 

It is important to compare the optimized focusing effect of the alter- 
nating-gradient system to that for the all-convergent lens system. We 
make the comparison on the alternating-gradient system formed by 
adding a divergent lens of equal focal length in between the lenses of an 
all-convergent lens system; then we have s = 2L. The optimized maxi- 
mum displacement due to input ray slope is 

r a 's = 2r„'L 

for the convergent lens system, and is 

3.33 r 'L 

for the alternating-gradient system. It is remarkable that the focusing 
effect of the alternating-gradient system is so nearly the same as that of 
the all-convergent lens system. In practice it may be advantageous to 
get the focusing action in a manner that inherently reverses itself period- 
ically. This analysis shows that such structures are nearly as effective as 
those wherein the focusing effect is always convergent. 

Fig. 3 shows that the focal length which is optimum with respect to 
the input ray slope (fc 2 ) is also an acceptable region with respect to input 
ray displacement (hi and fcs). 

VII. RAY PATHS 

One can get a useful physical feel for the wave propagation by tracing 
the rays in a few of the important cases. 

For the all-convergent lens system optimized according to (50), Fig. 
8 shows the ray paths for a zero-slope finite-displacement input ray and 
for a zero-displacement finite-slope input ray. Here 8 = 90° [see (35)] 
and a period is completed in 4 lenses. As proved in Appendix B, the 

* Note that k 9 = \/2, with s = 2/ in (36), but A- 9 cos (pS — fg) is always ±1 
for any p. (The angle ip 3 = 45°.) This is an example of the caution that must be 
exercised in regarding the A's as maximum values of the various terms in r„ , r m , 
r v , etc. 
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response to an arbitrary input ray can be obtained by a linear super- 
position of the responses shown in Fig. 8. 

For the alternating-gradient system the optimum according to (53) 
corresponds to an angle 6 in (15), (20), (26) and (30) of 76.4°, which 
makes the ray path periodic only at a very large number of lenses. How- 
ever, a very useful feel can be obtained from the ray paths for 6 = 90°, 
corresponding to L = y/2j and giving a value of k 2 only slightly larger 
than the minimum value (3.414 compared to 3.33). These ray plots are 
shown in Figs. 9 and 10 for the two types of input rays at the two pos- 
sible points in the alternating-gradient system. We note that input ray 
displacement causes the same maximum displacement in the response 
regardless of where it occurs. Input ray slope is much more serious when 



p=o 



pal 



p=2 



P = 3 



P=4 



P = 5 




Fig. 8 — Ray paths in the confocal all-convergent lens system, s = 2/. 

it occurs in front of a diverging lens than when it occurs in front of a 
converging lens. Again, the response to arbitrary input rays can be ob- 
tained by adding the plotted responses. 

One can gain a little feel for the stability comparison made previously 
by looking at Figs. 11 and 12. Even though the r term and the r</ term 
of (35) for the all-convergent lens system go to infinity individually 
when s = 4/, a suitable combination of input ray slope and displacement 
remains bounded and this is illustrated in Fig. 11.* Any reduction in 
focal length / causes instability, and any increase in / leaves the system 
completely stable. It is clear that adding a lens of any kind at the mid- 
point between lenses in Fig. 11 will not alter the propagation of that ray. 
In Fig. 12 we see that adding a divergent lens in between the converging 

* One can obtain these values of r, from (35) by a suitable limiting process. It 
is helpful to start with the alternative form of (35) : 



r P = r 



cos pS + .... . _ jk sin p&\ + r 's hi, sin pS. 
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SLOPE P' 




Fig. 9 — Ray paths for system of Fig. 1, L = y/2f. 

lenses will cause reductions in the focal length of the p = 1 converging 
lens to make the ray sent on to the p = 2 lens diverge even more; for 
increases in the focal length of the p = 1 lens, the divergent lens reduces 
the angle of the ray sent on to the p = 2 lens. Hence, it is plausible that 
the addition of the divergent lens between the convergent lenses does 
not alter the stability requirement on the focal lengths. 

Given the mathematically-derived condition that divergent lenses of 



TnL 
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Fig. 10 — Ray paths for system of Fig. 3, L = V'2f. 
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p = 4 




S -2 S 

RAY PATH FOR S = 4f 



pig |i _ R a y path for s = 4/ in all-convergent lens system. 

focal length / added to a chain of convergent lenses of focal length / do 
not change the stability criterion as described above, the physical argu- 
ment just outlined leads to the conclusion that weaker divergent lenses 
would also leave the stability criterion unchanged.* 



p=i 



P = 2 




Fig. 12 — Ray path for / near s/4 to illustrate effect of divergent lens. 



VIII. CONCLUSION 

Alternating-gradient focusing is surprisingly close to an all-convergent 
lens system in focusing ability, and may be preferred if practical matters 
such as structural features or cost favor the alternating-gradient system. 

* When reading this manuscript, Mr. J. P. Gordon commented that this con- 
clusion is in agreement with the work of Boyd and Kogelnik 3 which can be shown 
to yield the relation /a > /i - L/2 for the required focal length U of the diverg- 
ing lens in terms of the focal length /i of the converging lens and the spacing L. 
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APPENDIX A 



We examine here the maximum value that the term r 's fcio sin p8 of 
(35) can take as a function of lens number p when our objective is to 
minimize the term through appropriate choice of focal length. In the 
body of the article it has been shown that /, 1(l has a minimum at s = 2/. 
This corresponds to a value of 8 = tt/2 from (39) and it is evident that 
sin pd = sin p ir/2 is cither zero or unity for all integral values of p. 

In the more general case we want to know the value of 

/.m sin p8 = sin pfi/sin 8. (54) 

When it is recognized that p may take on all integral values greater than 
zero it follows that the maximum value of [(sin p8)/sm 8] as p varies 
can never be less than unity for any fixed 8. 

It is possible for [(sin p5)/sin 8) to have a maximum value which is 
smaller than k 10 = 1/sin 8. That is to say, sin p8 does not necessarily go 
through unity even though p ranges from to =c in integral steps. 

Referring to Fig. 13, the maximum value of sin p8 will be less than 
unity if 

S(8 + h) = *72 (55) 

or 

8 = Tr/(2q+l) (56) 

where q = 1, 2, 3, • • • .It also is true that sin p8 will have a maximum 
value less than unity for 

8 = r[*/(2q+ 1 )] (57) 

where r = 1, 2, 3, 4, • • • . 

The values of s/f corresponding to these values of 8 and the resultant 
values of maximum A'i sin p8 are given in Table I. Column 5 shows the 
ratio of A:i to the maximum of A- ]0 sin p8, and is a measure of the error 
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Fig. 13 — Diagram of sin pS yielding | sin pS \ < 1 for all p. 
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Table I 



1 

6 


2 

s/f 


kio 


4 

Max. Value of 

Aio sin pS 


5 

Column 3 + 

Column 4 


ir/2 

tt/3 

2tt/3 

7r/5 

2tt/5 
3tt/5 
4tt/5 
r/7 
2w/7 
3tt/7 
4tt/7 
5tt/7 
6ir/7 


2 
1 
3 

0.38 

1.38 

2. 62 

3.62 

0.194 

0.750 

1.554 

2.444 

3.226 

3.806 


1.0 

1.15 

1.15 

1.70 

1.05 

1.05 

1.70 

2.31 

1.28 

1.023 

1.023 

1.28 

2.31 


1.0 

1.0 

1.0 

1.62 

1.0 

1.0 

1.62 

2.255 

1.25 

1.0 

1.0 

1.25 

2.255 


1.0 

1.15 

1.15 

1.05 

1.05 

1.05 

1.05 

1.023 

1.023 

1.023 

1.023 

1.023 

1.023 



made in assuming sin p8 goes through unity. That ratio is 1/cos (5/2) 
where 8 is given by (56). All values for a given q in (57) result in the same 
error, but the various values of r indicate the values of 8 and s/f at which 
that error will appear. Fig. 14 summarizes the data of Table I; for p 
ranging up to infinity it is only at the discrete values of 8 given by (57) 
that the maximum of k i0 sin p8 differs from A' ]0 . 
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Fig. 14 — Maximum value of A- 10 sin pS vs s/f. 
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If p were finite and the ratio s/f was varied, the plot of Fig. 14 would 
presumably show finite-width dips of the same over-all depth as those 
plotted. 

APPENDIX B 

It is the purpose of this appendix to point out that the "thin lens" de- 
scription of light ray propagation leads to the following conclusion (see 
Fig. 15): 

The slope and displacement of the output ray of an arbitrary sequence 
of lenses for an input ray of slope r ' and displacement n is exactly the 
algebraic sum of the slopes and displacements of the output ray found 
(i) for an input ray of slope /■</ with zero displacement from the axis, 
and assuming all lens displacements d„ = 0, (ii) for an input ray of dis- 
placement Ti with zero slope and assuming all lens displacements d„ = 0, 
and (in) for an input ray of zero slope and zero displacement and as- 
suming one lens displacement at a time is nonzero, summing the ray 
output slopes and displacements thus found over all lens displacements. 

The proof is as follows: For the nth lens in a sequence of lenses (see 
Fig. 1): 



/ / _ ( r n - d„ \ 



(58) 



where /•„' is the slope of the ray immediately following the nth lens and 
r„ is the displacement at the nth lens. We note that the angular lens 
rotation tp n does not affect the ray propagation, an approximation which 
implies that 



or 



(r„ — d n ) cos ip„ ^ (/■„ — d H ) 



<Pn « 1. 



(59) 



LENS NO. ( NO. 2 

(FOCAL LENGTH =f,) f 2 



SLOPE =r' 




AXIS OF 
LENS SYSTEM 



- S 2 »+* S 3 * 

Fig. 15 — Light-ray path in an arbitrary lens system. 
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We may expand (58) to form 

, _ / _ (n - di) _ (r 2 - ja) _ (r„ - d n ) , QO j 

T n — To « f f ' 

./i h J" 

Similarly, the displacement /•„ at the rath lens is 

/•„ = ,-„_, + w-1' ( 61 ) 

which may he expanded to 

(r, - di) 



r„ = r„_i + s„_i <r ' - 



/i 
__ (r2 — <h) _ _ (r„_ x — d n -i) 

jl fn-l 



which is valid for n ^ 2. 

We may examine each term of (60) and (02) and find that 



(62) 



/ = A n r»' + B n n + "Z<* m ^ (63) 



r n = C»ro' + D n n + Eft»r (64) 

in which A„ , B„ , C n , D„ and the a,„ and /3„, are all independent of 
r a ', ri and the d m . 

We have thus proven the above-stated conclusion. 

REFERENCES 

1. Pierce, J. R., Theory and Design of Electron Beams, 2nd ed., D. Van Nostrand, 

Princeton, 1954. 

2. Courant, E. D., Livingston, M. S., and Snyder, H. S., The Strong-Focusing 

Synchrotron — A New High-Energy Accelerator, Phys. Rev., 88, Dec. 1, 
1952, pp. 1190-119(1. _, 

3 Boyd G D., and Kogelnik, H., Generalized Confocal Resonator Iheory, 
B.S.T.J., 41, July, 1962, pp. 1347-1369. 



